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Abstract. - Strong dynamical scaling violations exist in quenched two-dimensional systems 
with vector 0(3) order parameters. These systems support non-singular topologically stable 
configurations (slcyrmions). By tuning the stability of isolated skyrmions to expand or shrink, 
we find dramatic differences in the dynamical multiscaling spectrum of decaying moments 
of the topological charge density distribution and in particular in the decay of the energy-density 
e ~ (IpD- We present a simple two- length-scale model for the observed exponents in the case 
when isolated skyrmions expand. No such simple model is found when isolated skyrmions shrink. 



Phase-ordering systems quenched from disordered initial conditions into an ordered phase 
typically evolve correlations that dynamically scale — they are self-similar in time Dy- 
namical scaling provides a powerful framework for the analysis of phase-ordering exponents 
and correlations. To understand scaling, systems that violate scaling are important to study. 
Logarithmic scaling violations exist in the spherical limit of conserved vector systems 
Recently, strong scaling violations have been found in one-dimensional Q and two-dimensional 
(2D) iQ, 1^ systems supporting non-singular topological textures ("skyrmions"). 

In this paper we investigate the phase-ordering process of the 2D 0(3) vector model in 
more detail, paying particular attention to sub-dominant terms in the free energy. These terms 
control the stability of individual skyrmions towards shrinking or growing in size. Surprisingly, 
these sub-dominant terms dramatically change the exponent c/jn characterizing the decay of 
energy in the system. 

We restrict ourselves to dissipative noise-free "model- A" dynamics, dtrn = —SF/Srn, where 
I TO I — 1 is maintained as a constraint. Our effective coarse-grained free energy is 
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where the topological charge density p is 

p(r) — rh ■ {dxTTL x 9j,m)/47r, (2) 

and is locally conserved in continuum systems |^. The term tt9 J cPv in Eqn. (|^) is 
analogous, when 9 is positive, to the Skyrme term used to stabilize skyrmion configurations 
in high-energy models 0. At time t = 0, we randomly orient the unit vector m at each site. 
This initial condition corresponds to a quench from T = oo to T = 0. 

A 2D skyrmion (see, e.g., |^) is a field configuration rn{v) that wraps around the 0(3) 
order-parameter sphere exactly once. It is described by a spatially extended topological charge 
p(r) that integrates to ±1. Each skyrmion can be characterized by a size, defined, e.g., as the 
half- width of the typically bell-shaped spatial distribution of topological charge p(r). In a 2D 
0(3) phase-ordering system, which has a mixture of skyrmions and anti-skyrmions, we can 
define at least three natural length-scales: the average skyrmion size Lt ~ i'^^, the average 
separation of skyrmions Ljv i*^" , and a length characterizing positive and negative charge 
separation — for example the inverse interface density, Lc ^ t'^'^' , of p = contours. All three 
length-scales are found to have distinct growth exponents. 

The static (minimal energy) configurations of an infinite continuum system with a pure 
quadratic gradient energy (9 — 0) were studied by Belavin and Polyakov (BP) ||]. For any 
net number of skyrmions in the system, they found a class of minimal energy states that are 
degenerate with respect to both the scale and position of the individual skyrmions. Essentially, 
the skyrmions in the pure skyrmion (or pure anti-skyrmion) BP configurations do not interact, 
and they neither expand nor shrink. 

In generic systems, the energy density will also contain higher-order gradient terms, arising 
e.g. from lattice interactions. The dominant corrections at late times (and hence large length- 
scales) in quenched systems will come from fourth-order gradient terms, such as the Skyrme 
term in Eqn. (|^). Such terms will destabilize the individual skyrmions towards expanding or 
shrinking. By adding the Skyrme term in Eqn. (|l]) with a sufficiently large amplitude 9 > 9c, 
it is possible to ensure that all skyrmions expand (cl). In our simulations, we use a 9-point 
Laplacian with nearest and next-nearest neighbor interactions and an isotropic fourth-order 
term. Our leading interactions in Fourier space are proportional to (fc^ — fc^/24)|mj^p, resulting 
in a positive 9c- 

Previous work |^, ||, |lOj has considered systems where isolated skyrmions shrink and 
ultimately unwind through the lattice {9 < 9c), thus violating charge conservation. For 9 > 9c, 
topological charge density is eliminated solely through the mutual annihilation of regions of 
positive and negative p — a process conserving the topological charge. For 9 < 9c, skyrmion 
unwinding and skyrmion/anti-skyrmion annihilation are competing processes — a partially 
annihilated skyrmion cannot unwind by itself We re-analyze our earlier results j5|, and 
also study systems with 9 > 9c using the same numerical techniques. We use 9 = 0.012 
in dimensionless units with a square-lattice spacing of 0.01. Our results are qualitatively 
unchanged for 9 = 0.006, with identical exponents. This indicates that the value 9 = 0.012 is 
well above 9c- The results presented below were obtained in a system of size 512 x 512, and 
averaged over 10 (resp. 12) runs in the case 9 — (resp. 9 — 0.012). 



(^) Our Skyrme term in Eqn. ^ does not describe all possible fourth-order gradient terms. For 
general fourth-order terms, the energy of a spherically symmetrical skyrmion of size R is E{t)_^ 
1 -|- A{6 — 9c)R~^ , where ^ is a positive constant. The rate of change of the energy is dtE ~ B? 
Self-consistency gives dtR ~ —R~^, which leads (in our case) to R{t) cx [{6 — 6c){t — to)Y^'^ - Skyrmions 
are also destabilized towards shrinking, with a similar rate, in the continuum model when a soft-spin 
potential that allows variations in the spin magnitude is used 0. 
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To characterize the system in the absence of dynamical scaling of correlations, wc measure 
the asymptotic (late-time) decay exponents /3„ of the moments of the topological charge density 
distribution P{\p\,t): 

(|p(r,t)r)^t-'3", (3) 
where (• • •) denotes a spatial average. [We assume and observe power-law decay with time. 



The moment analysis in terms of /9„ is analogous to multifractal analysis of static fractal 1 12 
or turbulent [ p^ systems, only with the measuring length-scale replaced by the inverse time 
elapsed since the quench (see, e.g., Q). 

The results for the asymptotic evolution of the energy density e are dramatically different 
for 6 < 6c and 9 > 9c- At late times, e is dominated by the total number density of skyrmions 
and anti-skyrmions, so that e ~ {\p\) ^ ~ t~'^'pN ^nd /3i = 2(j>N. Whereas systems with 

skyrmions unstable towards shrinking {9 < 9c) have (pN = 0.32±0.01, systems with 9 > 9c have 
(pM = 0.23 ± 0.01 , as shown in Fig. 1. [Here we extract (t)N from the exchange energy, which 
dominates the total energy at late times (see Fig. 1). The values extracted directly from 
are (f)N — 0.31 ± 0.01 and (pN = 0.21 ± 0.01, respectively, and are consistent. Error bars shown 
are always statistical errors.] To our knowledge this is the first example of a phase-ordering 
system in which a sub-dominant interaction controls the growth laws. 

Higher moments of the topological charge density, Eqn. (|3|), are shown in Fig. 2. The /3„ 
are multiscaling — they are not simply proportional to n. This indicates that i^d/ol, i), shown 
in Fig. 3, is determined by more than one length scale. 

We can explain many of the observed features of P{\p\^t) using the profile of a single 
Belavin-Polyakov (BP) skyrmion in combination with the multiple length-scales L^j and 
Lc- We first discuss the case 9 > 9c, where (after initial transients due to large gradients) 
skyrmions grow and annihilate with anti-skyrmions, and there is no unwinding through the 
lattice. 

For 9 > 9c, Si scaling argument involving both Lt and Ln suffices to characterize all 
of the Pn- Picture the system divided into regions of size L^, each containing a skyrmion 
or anti-skyrmion of size Lt (see figures in [^, |^). Represent each skyrmion as an isolated 
BP configuration of scale Lt, so that p(r) — 7r~^Ly/(Ly + r"^)"^ [^ and (|p(r,t)l") oc 
Ln '^L'^^ J^"^ rdr [r'^ + L^)^^". For n > 1/2, the integral is dominated by the central region 
of radius r < Lt, and we obtain (|p(r, t)|") ^"ijv^. For n < 1/2, on the other hand, 

the integral is dominated by the large-r regions, giving (|p(r,t)|") ~ For n = 1/2, 

{\p{r, t)\^^^} ~ LtL^ \og{t). In addition, the mixture of skyrmions and anti-skyrmions in the 
quenched system requires the presence of contours on which p = (probed by the third length, 
Lc, which is discussed later). We observe P(|p|) = const for small which corresponds to a 
linear profile p{z) ^ z with distance z away from the p = contours. This constant regime for 
small p implies that moments with n < —I diverge, and the corresponding /3„ are not defined. 
Summarizing: 

r 2n(PT + 2{(j)N-(pT), n>l/2, 
'"~\nmN-2q}T), -l<n<l/2. ^> 

This piecewise-linear prediction agrees well with our data for 6* > Qj: in Fig. 2. The best fit is 
obtained for the values (f>N = 0.210 ± 0.005 and (f>T = 0.18 ± 0.02 (B). 

For 9 > 9c, we show our data for the distribution P{\p\,t) in Fig. 3 (lower set of curves). 
The P{\p\) = const regime is visible for small \p\, followed by a developing power-law regime 



(^) These exponents are distinct, as can be seen from the distinct linear regimes shown in the inset 
of Fig. 2. In light of this, there is covariance in the measured errors. 
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for intermediate \p\ and an exponential cutoff at large \p\. In Fig. 4, the data from Fig. 3 are 
rescaled to test the scaling ansatz t) — fi{t)G[\p\f2{t)], where /i and /2 are functions of 

time only. A good quality of collapse is achieved for large and intermediate a;, where x — \p\f2- 
However, significant scaling violations are seen at small x. 

The difitinct growth laws for L^r and Lt are sufficient to introduce scaling violations in 
i-*[|/9|, t\ (□). The approximate power-law behavior G[x] ^ x^^l"^ observed in the intermediate- a; 
range (see Fig. 4) is consistent with a BP-like profile of the skyrmions. For a skyrmion of size 
Lt, p(y) = 7r^^L|,/(L|, -I- r^)^ implies the topological charge distribution of a single skyrmion 

Plt{p) — ^LtP~^^'^ for p < pmax — TT~^ L^^ . However, the tail of the BP profile is expected 
to be cut off at r > Ln due to the presence of neighboring skyrmions and antiskyrmions; this 
corresponds to p values p > Pc ^ ^20r-40jv_ sj^ce Pmax/Pc ~ i4(0„-</,T)^ ^he width of the 
G{x\ ~ a;""^/^ regime is predicted to increase with t — hence violating the scaling ansatz. Our 
data is consistent with this, and the lack of scaling of P{\p\,t) is directly confirmed by the 
observed multiscaling character of (in- 

For 9 < Oc^ the simple picture presented above does not apply. The instability of skyrmions 
towards unwinding in this case leads to the presence of small skyrmions that significantly 
change P(lyol,t) and hence /3„. Furthermore, we show below that unwinding skyrmions cannot 
be treated independently of annihilation. To proceed, we consider the distribution density 
H{R,t) of skyrmion sizes R at time t (□). 

We see in Fig. 2 that /3„ saturate asymptotically at large n, for 6 = 0. This indicates the 
presence of small unwinding skyrmions, which dominate the moments of P(|/9|,i) for large 
enough n. Let us first assume that small skyrmions (with sizes R Lj-) do not participate in 
annihilation processes and unwind independently. In that case, the size distribution function 
iJ(i?, t) satisfies a continuity equation dtH + duJ = 0, where the "size current" is (see footnote 
'^^^) J = HdtR ~ —H/R^. To maintain N , where N is the total number density of skyrmions 
and antiskyrmions, we must have a constant flux J{R) = J{C) ^ N (where ^ is the lattice 
scale), OT H —NR"^ for small R. These small skyrmions would dominate the high /?„, so 
that (IpI") - /^^" dRH{R)R'^-^'' - Ar^6-2« _ ^-(i+20„) „ > 3_ 

For 9 = 0, the measured /3„ do approach 1 + 2(j)N. This indicates that unwinding is a 
signiflcant relaxation process at late times. However, they are not constant for large n but 
instead fit the empirical form 

/3„ = 1 + 2(bN - n-P (5) 

within error bars for all n > 1 (the unit prefactor of the power-law term ensures that (3i = 2(1) n)- 
In Fig. 2, we plot the best fit, with (jyjy = 0.31 and p = 0.78 {(pN — 1/3 and p = 2/3 provide 
a comparable fit). Our data for P{\p\,t) arc shown in Fig. 3 for = (upper set of curves). 
Similarly to the 6 > 6^ case, we see a P{\p\) = const regime at low \p\. As seen in Fig. 4, 
P{\p\,t) does not scale well, even for large x. An asymptotic power-law tail P{\p\,t) ~ 
would lead to /3„ saturating for n > 3. The asymptotic approach towards that power-law 
(shown in Fig. 3) is another manifestation of /3„ asymptotically approaching 1 -f 2(/>jv at large 
n. From the previous paragraph, the existence of the non-zero correction term n~P for n > 3 
is inconsistent with having a sub-population of independently unwinding BP skyrmions. 

Two equivalent mechanisms can describe this breakdown of the picture of small unwinding 
BP skyrmions for 9 < 9c- Either the profile of individual shrinking skyrmions is asymmetric 
and hence not BP like, or the partial annihilation of skyrmions with anti-skyrmions changes 



(^) Additional scaling violations arise at small x arise due to a third length-scale Lc- 
To obtain Eqn. (g) we eflFectively used H{R,t) = 5[R- Lrit)] IL\. 
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the unwinding process (13). In comparison, high-order gradient terms on their own do not 
significantly affect the structure of an unwinding skyrmion until it approaches the lattice 
scale. 

On the other hand, does unwinding affect the structure and evolution of annihilating 
skyrmions and anti-skyrmions? To address this question, we need to probe the /an7e-scale 
spatial distribution of the skyrmions and anti-skyrmions. We can do this by measuring the 
interface line-density of p = 0. We measure this density, L^^ t~'^'^ , by counting the nearest 
neighbor lattice-sites between which p changes sign, and find (j)c = 0.41 ± 0.01 for both < 6^ 
and 6 > Of,. This exponent is significantly greater than (f)T or and is another demonstration 
of the strong scaling violation seen in our system. 

In addition to Lc, we have measured two other quantities that probe the large-scale 
structure of the system. The \ow-p value of the topological charge density P(0,t) ~ ^o.75±o.05 
probes the charge distribution near the p — contours, while the half-width Lp ^ ^o.22±o.03 
of p — p correlations provides a measure of the typical skyrmion size. The agreement of the 
given values of the exponents of Lc, -P(0, t), and Lp in the 9 > 9c and 9 < 9c cases is striking 
— especially considering the dramatic differences in the spectrum. 

We conclude simply: annihilation affects unwinding, but not vice versa. 

In summary, we have studied phase-ordering dynamics in 2D 0(3) systems supporting 
skyrmions. By tuning the strength 9 of the Skyrme term, we made isolated skyrmions unstable 
to shrink and unwind {9 < 9c) or to expand [9 > 9c). We found dramatic changes in the 
multiscaling spectrum /?„ of moments of the topological charge density between the two cases. 
Having expanding skyrmions leads to a piecewise linear /5„ spectrum explained by taking into 
account that the average skyrmion size Lxit) and the average skyrmion separation Liq(t) scale 
differently in time. Having shrinking skyrmions leads to a curved /3„ spectrum remarkably well 
fit by Eqn. (|). 

ADR acknowledges support by EPSRC grant GR/J78044; MZ by NSF grants DMR-89- 
20538, DMR-94-24511, and DMR-95-07366. We thank Shivaji Sondhi for discussions. 
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Fig. 1. - The exchange energy [first term of Eqn. (|])] for Skyrme-term amplitudes 9 — 0.012 
(upper curve) and = (middle curve). The lower curve shows the Skyrme energy [second 
term of Eqn. (|^)] in the case 9 = 0.012. Even though the exchange energy dominates the total 
energy at late times, the decay exponent of the exchange energy depends on the Skyrme-term 
amplitude 9. The straight dot-dashed lines have slopes —0.46 (upper line), —0.64 (middle), 
and —0.78 (lower). 
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Fig. 2. - Growth exponents /?„ for the n^^ moment of the topological charge density, Eqn. (||). 
Shown are the /3„ spectra for a system with Skyrme-term amplitude = < 9c and with 
9 = 0.012 > dc- Error bars do not exceed the size of the symbols. The full line shows the 
two-length-scale prediction of Eqn. (Q) with 0t = 0.187 and (/)n = 0.210. The dotted-dashed 
line shows the phenomenological fit to the /?„ spectrum in the 9 = case, given by Eqn. (|^) 
with = 0.31 and p = 0.78. Note that this fit is incompatible with the data for n < 1. 
The inset highlights more clearly the discontinuity a,t n = 1/2 in the slope of the /3„ curve 
for 9 = 0.012, where the dashed lines are linear extrapolations of the regions above and below 
n = 1/2. 
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Fig. 3. - The distribution P{\p\,t) of topological charge density at indicated times t. The 
evolution of the small p regime is well described by P(0, t) ~ f0.75±o.05 ^^^^^ 

cases. The 

vertical axis of the 9 = 0.012 curves is offset by three decades for clarity. 
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Fig. 4. - Re-scaled distributions G{x, t) of topological charge density at indicated times t. The 
vertical axis of the 9 = 0.012 curves is offset by one decade for clarity. The horizontal and 

vertical axes of the original P(|p|,f) curves have been rescaled so that all curves collapse at 
X = 1. The straight lines drawn in the figure have slopes —3/2 (dashed line) and —4 (solid 
line). 



